BNL-NT-02/5 



Classical Chromo-Dynamics of 
Relativistic Heavy Ion Collisions * 

Dmitri Kharzeev 

Physics Department, Brookhaven National Laboratory, Upton, NY '11973-5000 
March 15, 2002 

Abstract. Relativistic heavy ion collisions produce thousands of particles, and it is 
sometimes difficult to believe that these processes allow for a theoretical description 
directly in terms of the underlying theory - QCD. However once the parton densities 
are sufficiently large, an essential simplification occurs - the dynamics becomes semi- 
classical. As a result, a simple ab initio approach to the nucleus-nucleus collision 
dynamics may be justified. In these lectures, we describe the application of these 
ideas to the description of multi-particle production in relativistic heavy ion colli- 
sions. We also discuss the role of semi-classical fields in the QCD vacuum in hadron 
interactions at low and high energies. 
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1. What is Chromo Dynamics? 

Strong interaction is, indeed, the strongest force of Nature. It is re- 
sponsible for over 80% of the baryon masses, and thus for most of the 
mass of everything on Earth and in the Universe. Strong interactions 
bind nucleons in nuclei, which, being then bound into molecules by 
much weaker electro-magnetic forces, give rise to the variety of the 
physical World. Quantum Chromo-Dynamics is the theory of strong 
interactions, and its practical importance is thus undeniable. But QCD 
is more than a useful tool - it is a consistent and very rich field theory, 
which continues to serve as a stimulus for, and testing ground of, many 
exciting ideas and new methods in theoretical physics. 

These lectures will deal with QCD of strong color fields, which can 
be explored in relativistic heavy ion collisions. (See the lectures by E. 
Iancu, A. Leonidov, L. McLerran (Iancu, Leonidov, McLerran, 2002), 
A.H. Mueller (Mueller, 2001), and R. Venugopalan in this volume for 
complementary presentation of the subject and more details.) 

* Lectures given at Cargese Summer School on QCD Perspectives on Hot and 
Dense Matter, Cargese, France, 6-18 Aug 2001. 
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1.1. QCD: THE Lagrangean 

So what is QCD? From the early days of the accelerator experiments 
it has become clear that the number of hadronic resonances is very 
large, suggesting that all hadrons may be classified in terms of a smaller 
number of (more) fundamental constituents. A convenient classification 
was offered by the quark model, but QCD was not born until the 
hypothetical existence of quarks was not supplemented by the prin- 
ciple of local gauge invariance, previously established as the basis of 
electromagnetism. The resulting Lagrangian has the form 

L = -]G%G% + J2 m-lvD* ~ m f )q a f, (1) 
/ 

the sum is over different colors a and quark flavors /; the covariant 
derivative is = — igA^t®, where t a is the generator of the 
color group SU(3), is the gauge (gluon) field and g is the coupling 
constant. The gluon field strength tensor is given by 

G% = d^A a v - d v Al + <?r bc 4^, (2) 
where f abc is the structure constant of SU(3): [t a ,t b ] = if abc t c . 

1.2. Asymptotic freedom 

Due to the quantum effects of vacuum polarization, the charge in field 
theory can vary with the distance. In electrodynamics, summation 
of the electron-positron loops in the photon propagator leads to the 
following expression for the effective charge, valid at r » vq: 

This formula clearly exhibits the "zero charge" problem (Landau and 
Pomeranchuk, 1955) of QED: in the local limit tq — ► the effective 
charge vanishes at any finite distance away from the bare charge due 
to the screening. Fortunately, because of the smallness of the physical 
coupling, this apparent inconsistency of the theory manifests itself only 
at very short distances ~ exp{— 37r/[2a em ]}, a em ~ 1/137. Such short 
distances are (and probably will always remain) beyond the reach of 
experiments, and one can safely use QED as a truly effective theory. 

As it has been established long time ago (Gross and Wilczek; Politzer, 
1973), QCD is drastically different from electrodynamics in possessing 
the remarkable property of "asymptotic freedom" - due to the fact 
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that gluons carry color, the behavior of the effective charge a s = g 2 /4ir 
changes from the familiar from QED screening to anti-screening: 



as long as the number of flavors does not exceed 16 (N c = 3), the 
anti-screening originating from gluon loops overcomes the screening 
due to quark-antiquark pairs, and the theory, unlike electrodynamics, 
is weakly coupled at short distances: a s (r) — ► when r — > 0. 

1.3. Chiral symmetry 

In the limit of massless quarks, QCD Lagrangian (1) possesses an ad- 
ditional symmetry Ul{N^) x U R (Nf) with respect to the independent 
transformation of left- and right-handed quark fields qL, R = ^±75)5: 



this means that left- and right-handed quarks are not correlated. Even 
a brief look into the Particle Data tables, or simply in the mirror, can 
convince anyone that there is no symmetry between left and right in 
the physical World. One thus has to assume that the symmetry (5) 
is spontaneously broken in the vacuum. The flavor composition of the 
existing eight Goldstone bosons (3 pions, 4 kaons, and the 77) suggests 
that the U A (1) part of U L (3) x U R (3) = SU L (3) x SU R (3) x U v (l) x 
UaO-) does not exist. This constitutes the famous "[/^(l) problem". 

1.4. The origin of mass 

There is yet another problem with the chiral limit in QCD. Indeed, as 
the quark masses are put to zero, the Lagrangian (1) does not contain 
a single dimensionful scale - the only parameters are pure numbers 
N c and Nf. The theory is thus apparently invariant with respect to 
scale transformations, and the corresponding scale current is conserved: 
du.Su = 0. However, the absence of a mass scale would imply that all 
physical states in the theory should be massless! 

1.5. Quantum anomalies and classical solutions 

Both apparent problems - the missing Ua(1) symmetry and the origin 
of hadron masses - are related to quantum anomalies. Once the cou- 
pling to gluons is included, both flavor singlet axial current and the 
scale current cease to be conserved; their divergences become propor- 
tional to the a s G" G° and gluon operators, correspond- 
ingly. This fact by itself would not have dramatic consequences if the 



a s (r) ~ 



{11NJ2- N f )\n(r / r y 



(4) 



QL V L q L ; q R -► V R q R ; V L , V R € U(N f ); 



(5) 
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gluonic vacuum were "empty", with G° = 0. However, it appears 
that due to non-trivial topology of the SU(3) gauge group, QCD 
equations of motion allow classical solutions even in the absence of 
external color source, i.e. in the vacuum. The well-known example of 
a classical solution is the instanton, corresponding to the mapping of 
a three-dimensional sphere S 3 into the SU(2) subgroup of SU(3); its 
existence was shown to solve the ?7a(1) problem. 

1.6. Confinement 

The list of the problems facing us in the study of QCD would not be 
complete without the most important problem of all - why are the 
colored quarks and gluons excluded from the physical spectrum of the 
theory? Since confinement does not appear in perturbative treatment of 
the theory, the solution of this problem, again, must lie in the properties 
of the QCD vacuum. 

1.7. Understanding the Vacuum 

As was repeatedly stated above, the most important problem facing us 
in the study of all aspects of QCD is understanding the structure of 
the vacuum, which, in a manner of saying, does not at all behave as 
an empty space, but as a physical entity with a complicated structure. 
As such, the vacuum can be excited, altered and modified in physical 
processes (Lee and Wick, 1974). 



2. Strong interactions at short and large distances 

In this lecture we will investigate the influence of QCD vacuum on 
hadron interactions at short and large distances. To make the problem 
treatable, we will limit ourselves to heavy quarkonia. In this lecture 
I will describe two recent results - one on the scattering of heavy 
quarkonia at very low energies, another on high-energy scattering. The 
common idea behind these two examples is to explore the influence 
of the QCD vacuum on hadron interactions. The presentation will be 
schematic, and I refer the interested reader to the original papers (Fujii 
and Kharzeev, 1999) and (Kharzeev and Levin, 2000) for details. 
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2.1. The long-range forces of QCD 



2.1.1. Perturbation theory 

Let us begin with a somewhat academic problem - the scattering of 
two heavy quarkonium states at very low energies. The Wilson operator 
product expansion allows one to write down the scattering amplitude 
(in the Born approximation) of two small color dipoles in the following 
form(Bhanot and Peskin, 1977): 

V(R) = -ij dt(0\T (5>Oi(0)) (^cjOjix^j |0>, (6) 

where x = (t,R), Oi(x) is the set of local gauge-invariant operators 
expressible in terms of gluon fields, and q are the coefficients which 
reflect the structure of the color dipole. At small (compared to the 
binding energy of the dipole) energies, the leading operator in (6) is the 
square of the chromo-electric field (l/2)g 2 E 2 (Voloshin, 1978; Bhanot 
and Peskin, 1977). Keeping only this leading operator, we can rewrite 
(6) in a simple form 

V(R) = -*(d 2 ^) 2 j dt(0\T ^ 2 E 2 (0) ±g 2 E\t,R)\0), (7) 

where d~2 is the corresponding Wilson coefficient defined by 

where we have explicitly factored out the dependence on the quarko- 
nium Bohr radius ao and the Rydberg energy eo; N is the number of 
colors, and \(f>) is the quarkonium wave function, which is Coulomb 
in the heavy quark limit 1 . In physical terms, the structure of (7) is 
transparent: it describes elastic scattering of two dipoles which act on 
each other by chromo-electric dipole fields; color neutrality permits 
only the square of dipole interaction. It is convenient to express <? 2 E 2 
in terms of the gluon field strength tensor (Novikov and Shifman, 1981): 

g 2 E 2 = - l -g*G aP G a P +g 2 (-G 0a G% + \g^G a(3 G af3 ) = 



\ ' <f<C 00 



b 



1 The Wilson coefficients c?2, evaluated in the large N limit, are available for S 
(Bhanot and Peskin, 1977) and P (Kharzeev, 1996) quarkonium states. 
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where 

f 



9„ = ^G afia Gl R = -^G af3a G a nR . (10) 



Note that as a consequence of scale anomaly, 6>^ is the trace of the 
energy-momentum tensor of QCD in the chiral limit of vanishing light 
quark masses. 

Let us now introduce the spectral representation for the correlator 
of the trace of energy-momentum tensor: 

<O|T0£(O)0£(s)|O> = J da 2 Pe (a 2 )A F (x;a 2 ), (11) 

where pe(a 2 ) is the spectral density and Ap(x;a 2 ) is the Feynman 
propagator of a scalar field. Using the representation (11) in (7), we 
get 

<*<*> - -(*f) 2 (?)7^>^r" fi (12) 

The potential (12) is simply a superposition of Yukawa potentials cor- 
responding to the exchange of scalar quanta of mass a. 

Our analysis so far has been completely general; the dynamics enters 
through the spectral density. In perturbation theory, for SU(N), one 
has 

pt, 2n ( bg 2 V N 2 -l 4 

Substituting (13) into (12) and performing the integration over invari- 
ant mass a 2 , we get, for N = 3 

The oc R~ 7 dependence of the potential (14) is a classical result 
known from atomic physics (Casimir and Polder, 1948); as is apparent 
in our derivation (note the time integration in (12)), the extra as 
compared to the Van der Waals potential oc i?~ 6 is the consequence 
of the fact that the dipoles we consider fluctuate in time, and the 
characteristic fluctuation time r ~ e$ 1 , is small compared to the spatial 
separation of the "onia" : r <C R. 

Let us note finally that the second term in (9) gives the contribution 
of the same order in g; this contribution is due to the tensor 2 ++ state 
of two gluons and can be evaluated in a completely analogous way. 
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(a) (b) 

Figure 1. Contributions to the scattering amplitude from (a) two gluon exchange 
and (b) correlated two pion exchange. 

Adding this contribution to (14), changes the factor of 15 in (14) to 23, 
and we reproduce the result of ref. (Bhanot and Peskin, 1977), which 
shows the equivalence of the spectral representation method used here 
and the functional method of ref. (Bhanot and Peskin, 1977). 

2.1.2. Beyond the perturbation theory: scale anomaly and the role of 
pions 

At large distances, the perturbative description breaks down, because, 
as can be clearly seen from (12), the potential becomes determined by 
the spectral density at small q 2 , where the transverse momenta of the 
gluons become small. At small invariant masses, we have therefore to 
saturate the physical spectral density by the lightest state allowed in 
the scalar channel - two pions. Since, according to (10), 9% is gluonic 
operator, this requires the knowledge of the coupling of gluons to pi- 
ons. This looks like a hopeless non-perturbative problem, but it can 
nevertheless be rigorously solved, as it was shown in ref. (Voloshin and 
Zakharov, 1980) (see also (Novikov and Shifman, 1981)). The idea is 
the following: at small pion momenta, the energy-momentum tensor 
can be accurately computed using the low-energy chiral Lagrangian: 

= _Q^a d n 7T a + 2m\'K a 'K a + ■■■ (15) 

Using this expression, in the chiral limit of vanishing pion mass one 
gets an elegant result (Voloshin and Zakharov, 1980) 

(Ol^G^GyiV)^. (16) 

Now that we know the coupling of gluons to the two pion state, the 
pion-pair contribution to the spectral density can be easily computed 
by performing the simple phase space integration, with the result 

PVtf) = (17) 
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which leads to the following long-distance potential (Fujii and Kharzeev, 
1999): 

(18) 

Note that, unlike the perturbative result which is manifestly ~ g 4 , the 
amplitude (18) is ~ g - this "anomalously" strong interaction is the 
consequence of scale anomaly 2 . 

While the shape of the potential in general depends on the spectral 
density, which is fixed theoretically only at relatively small invariant 
mass, the overall strength of the non-perturbative interactions is fixed 
by low energy theorems and is determined by the energy density of 
QCD vacuum. Indeed, in the heavy quark limit, one can derive the 
following sum rule (Fujii and Kharzeev, 1999): 

J d 3 R(v(R)-VP\R)) = (J 2 ^) 2 (^) 2 16|e mc |, (19) 

which expresses the overall strength of the interaction between two 
color dipoles in terms of the energy density of the non-perturbative 
QCD vacuum. 

2.1.3. Does a s ever get large? 

Asymptotic freedom ensures the applicability of QCD perturbation 
theory to the description of processes accompanied by high momentum 
transfer Q. However, as Q decreases, the strong coupling a s (Q) grows, 
and the convergence of perturbative series is lost. How large can a s get? 
The analyzes of many observables suggest that the QCD coupling may 
be "frozen" in the infrared region at the value {a s ) IR ~ 0.5 (see (Dok- 
shitzer, 1998) and references therein). Gribov's program (Gribov, 1999) 
relates the freezing of the coupling constant to the existence of massless 
quarks, which leads to the "decay" of the vacuum at large distances 
similar to the way it happens in QED in the presence of "supercritical" 
charge Z > 1/a. One may try to infer the information about the 
behavior of the coupling constant at large distances by performing the 
matching of the fundamental theory onto the effective chiral Lagrangian 
at a scale Q ~ 47r/„- ~ 1 GeV, at which the ranges of validity of 
perturbative and chiral descriptions meet (Fujii and Kharzeev, 1999). 
It is easy to see that in the chiral limit the matching of the potentials 

2 Of course, in the heavy quark limit the amplitude (18) will nevertheless vanish, 
since ao — * and eo — > oo. 
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(18) and (14) yields 3 the coupling constant which freezes at the value 



/ \ 6y/2 7T 

{a s ) IR - niVc _ 

numerically, for QCD with N c = 3 and Nf = 2 one finds (a s )m ~ 

0. 56. Note that this expression has an expected N c dependence in the 
topological expansion limit of N c — > oo, Nf/N c = const. 

Since the trace of the energy momentum tensor in general relativity 
is linked to the curvature of space-time, the matching procedure leading 
to Eq. (20) has an interesting geometrical interpretation: it corresponds 
to the matching, at a relatively large distance, of curved space-time of 
the fundamental QCD with the flat space-time of the chiral theory. 

2.2. High-energy scattering: scale anomaly and the "soft" 

POMERON 

In a 1972 article entitled "Zero pion mass limit in interaction at very 
high energies" (Anselm and Gribov, 1972), A. A. Anselm and V.N. 
Gribov posed an interesting question: what is the total cross section 
of hadron scattering in the chiral limit of m n — > 0? On one hand, as 
everyone believes since the pioneering work of H. Yukawa, the range 
of strong interactions is determined by the mass of the lightest meson, 

1. e. is proportional to ~ rn ^ K 1 ■ The total cross sections may then be 
expected to scale as ~ "V 2 j and would tend to infinity as m w — > 0. 
On the other hand, soft-pion theorems, which proved to be very useful 
in understanding low-energy hadronic phenomena, state that hadronic 
amplitudes do not possess singularities in the limit m n — > 0, and one 
expects that the theory must remain self-consistent in the limit of the 
vanishing pion mass. At first glance, the advent of QCD has not made 
this problem any easier; on the contrary, the presence of massless gluons 
in the theory apparently introduces another long-range interaction. 
Here, we will try to address this problem considering the scattering of 
small color dipoles. 

Again, perturbation theory provides a natural starting point. In the 
framework of perturbative QCD, a systematic approach to high en- 
ergy scattering was developed by Balitsky, Fadin, Kuraev and Lipatov 
(Balitsky, Fadin, Kuraev and Lipatov, 1977), who demonstrated that 
the "leading log" terms in the scattering amplitude of type (g 2 ln s) n 
(where g is the strong coupling) can be re-summed, giving rise to the 
so-called "hard" Pomeron. Diagramatically, BFKL equation describes 

3 The matching procedure of course can be performed directly for the correlation 
function of the energy-momentum tensor. 



Nf 



(20) 
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the t— channel exchange of "gluonic ladder" - a concept familiar from 
the old-fashioned multi-peripheral model. 

It has been found, however, that at sufficiently high energies the per- 
turbative description breaks down (Mueller, 1997), (Dokshitzer, 1998). 
The physical reason for this is easy to understand: the higher the energy, 
the larger impact parameters contribute to the scattering, and at large 
transverse distances the perturbation theory inevitably fails, since the 
virtualities of partons in the ladder diffuse to small values. At this 
point, the following questions arise: Does this mean that the problem 
becomes untreatable? Does the same difficulty appear at large distances 
in low-energy scattering? And, finally, what (if any) is the role played 
by pions? 

The starting point of the approach proposed in (Kharzeev and Levin, 
2000) is the following: among the higher order, 0(a 2 s ) (as = g 2 /Air) 
, corrections to the BFKL kernel one can isolate a particular class of 
diagrams which include the propagation of two gluons in the scalar color 
singlet channel J PC = ++ . We then show that, as a consequence of 
scale anomaly, these, apparently 0(a s ), contributions become the dom- 
inant ones, 0(a s ). This is similar to our previous discussion in Sect. 2.1, 
where the interaction potential, proportional to a s , at large distances 
turned into a "chiral" potential ~ a s due to the scale anomaly. 

One way of understanding the disappearance of the coupling con- 
stant in the spectral density of the g 2 G 2 operator is to assume that 
the non-perturbative QCD vacuum is dominated by the semi-classical 
fluctuations of the gluon field. Since the strength of the classical gluon 
field is inversely proportional to the coupling, G ~ l/g, the quark 
zero modes, and the spectral density of their pionic excitations, appear 
independent of the coupling constant. 

The explicit calculation using the methods of (Gribov, Levin and 
Ryskin, 1983) yields the power-like behavior of the total cross section: 



oo 

^ \ " BORN ^ A foil 

Otot = 2^ a n = CP S , (21) 

n=0 



where o BOHN is the cross section due to two gluon exchange, and 
the non-perturbative contribution to the intercept A is (Kharzeev and 
Levin, 2000) 

A - t x if)' * 3^ / 1£ i^ M2 > ~ ^ CD ^))- 

(22) 
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Using the chiral formula (17) for pf ys for M 2 < M 2 , we obtain the 
following result (Kharzeev and Levin, 2000): 

The precise value of the matching scale Mq as extracted from the low- 
energy theorems depends somewhat on detailed form of the spectral 
density, and can vary within the range of Mq = 4^-6 GeV 2 . Fortunately, 
the dependence of Eq. (23) on Mq is only logarithmic, and varying it 
in this range leads to 

A = 0.08 -T- 0.1, (24) 

in agreement with the phenomenological intercept of the "soft" Pomeron, 
A ~ 0.08. 

At present, the language used in the description of hadron interac- 
tions at low and high energies is very different. Yet, as the two examples 
discussed above imply, both limits may appear to be determined by the 
same fundamental object - the QCD vacuum. 



3. QCD in the classical regime 

Most of the applications of QCD so far have been limited to the short 
distance regime of high momentum transfer, where the theory becomes 
weakly coupled and can be linearized. While this is the only domain 
where our theoretical tools based on perturbation theory are adequate, 
this is also the domain in which the beautiful non-linear structure 
of QCD does not yet reveal itself fully. On the other hand, as soon 
as we decrease the momentum transfer in a process, the dynamics 
rapidly becomes non-linear, but our understanding is hindered by the 
large coupling. Being perplexed by this problem, one is tempted to 
dream about an environment in which the coupling is weak, allowing 
a systematic theoretical treatment, but the fields are strong, revealing 
the full non-linear nature of QCD. I am going to argue now that this 
environment can be created on Earth with the help of relativistic heavy 
ion colliders. Relativistic heavy ion collisions allow to probe QCD in the 
non-linear regime of high parton density and high color field strength. 

It has been conjectured long time ago that the dynamics of QCD in 
the high density domain may become qualitatively different: in parton 
language, this is best described in terms of parton saturation (Gribov, 
Levin and Ryskin, 1983; Mueller and Qiu, 1986; Blaizot and Mueller, 
1987), and in the language of color fields - in terms of the classi- 
cal Chromo-Dynamics (McLerran and Venugopalan, 1994); see the 
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ln(l/x) 



Figure 2. The place of relativistic heavy ion physics in the study of QCD; the vertical 
axis is the product of atomic numbers of projectile and target, and the horizontal 
axes are the momentum transfer Q 2 and rapidity y = ln(l/a;) (x is the Bjorken 
scaling variable). 



lectures (Iancu, Leonidov, McLerran, 2002) and (Mueller, 2001) and ref- 
erences therein. In this high density regime, the transition amplitudes 
are dominated not by quantum fluctuations, but by the configura- 
tions of classical field containing large, ~ l/a s , numbers of gluons. 
One thus uncovers new non-linear features of QCD, which cannot be 
investigated in the more traditional applications based on the pertur- 
bative approach. The classical color fields in the initial nuclei (the 
"color glass condensate" (Iancu, Leonidov, McLerran, 2002)) can be 
thought of as either perturbatively generated, or as being a topolog- 
ically non-trivial superposition of the Weizsacker- Williams radiation 
and the quasi-classical vacuum fields (Kharzeev, Kovchegov and Levin, 
2001; Nowak, Shuryak and Zahed, 2001; Kharzeev, Kovchegov and 
Levin, 2002). 

3.1. Geometrical arguments 

Let us consider an external probe J interacting with the nuclear target 
of atomic number A. At small values of Bjorken x, by uncertainty 
principle the interaction develops over large longitudinal distances z ~ 
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1/mx, where m is the nucleon mass. As soon as z becomes larger 
than the nuclear diameter, the probe cannot distinguish between the 
nucleons located on the front and back edges of the nucleus, and all 
partons within the transverse area ~ 1/Q 2 determined by the momen- 
tum transfer Q participate in the interaction coherently. The density 
of partons in the transverse plane is given by 



PA 



xGa{x, Q 2 



A 1 / 3 , 



(25) 



where we have assumed that the nuclear gluon distribution scales with 
the number of nucleons A. The probe interacts with partons with cross 
section a ~ a s /Q 2 ; therefore, depending on the magnitude of momen- 
tum transfer Q, atomic number A, and the value of Bjorken x, one may 
encounter two regimes: 

— apA C 1- this is a familiar "dilute" regime of incoherent inter- 
actions, which is well described by the methods of perturbative 
QCD; 

— apA 3> 1 - in this regime, we deal with a dense parton sys- 
tem. Not only do the "leading twist" expressions become inade- 
quate, but also the expansion in higher twists, i.e. in multi-parton 
correlations, breaks down here. 



~ 1/Q 




Figure 3. Hard probe interacting with the nuclear target resolves the transverse 
distance ~ \j\fQ (Q 2 is the square of the momentum transfer) and, in the target 
rest frame, the longitudinal distance ~ 1/ (mi) (m is the nucleon mass and x - 
Bjorken variable). 



The border between the two regimes can be found from the condition 
a pa — 1; it determines the critical value of the momentum transfer 
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("saturation scale" (Gribov, Levin and Ryskin, 1983)) at which the 
parton system becomes to look dense to the probe 4 : 

2 xG A (x,Q 2 ) 
Q ° ~ °* ,R\ ■ (26) 

In this regime, the number of gluons from (26) is given by 

xG A {x,(fi)~-^^Q 2 a R 2 A , (27) 

where Q 2 R 2 A ~ A. One can see that the number of gluons is propor- 
tional to the inverse of a s (Q 2 ), and becomes large in the weak coupling 
regime. In this regime, as we shall now discuss, the dynamics is likely 
to become essentially classical. 

3.2. Saturation as the classical limit of QCD 

Indeed, the condition (26) can be derived in the following, rather gen- 
eral, way. As a first step, let us re-scale the gluon fields in the La- 
grangian (1) as follows: A® — > A® = gA a . In terms of new fields, 

G^ v = gG a ^ v = d^A a u - d v Al + f^A^Af,, and the dependence of the 
action corresponding to the Lagrangian (1) on the coupling constant is 
given by 

S ~ / \ d ' x - ( 28 ) 

Let us now consider a classical configuration of gluon fields; by defini- 
tion, Gy V in such a configuration does not depend on the coupling, and 
the action is large, S 3> h. The number of quanta in such a configuration 
is then 

where we re- wrote (28) as a product of four-dimensional action density 
P4 and the four-dimensional volume V4. 

Note that since (29) depends only on the product of the Planck 
constant h and the coupling g 2 , the classical limit h — > is indistin- 
guishable from the weak coupling limit g 2 — ► 0. The weak coupling 
limit of small g 2 = 47tq s therefore corresponds to the semi-classical 
regime. 

The effects of non-linear interactions among the gluons become im- 
portant when d^A^ ~ A 2 ^ (this condition can be made explicitly gauge 

4 Note that since xGa{x,Q 2 s ) ~ A 1 ^ 3 , which is the length of the target, this 
expression in the target rest frame can also be understood as describing a broadening 
of the transverse momentum resulting from the multiple re-scattering of the probe. 
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invariant if we derive it from the expansion of a correlation function 
of gauge-invariant gluon operators, e.g., G 2 ). In momentum space, this 
equality corresponds to 

Q 2 ~A 2 ~ (G 2 ) 1/2 = Vpi; (30) 

Q s is the typical value of the gluon momentum below which the inter- 
actions become essentially non-linear. 

Consider now a nucleus A boosted to a high momentum. By uncer- 
tainty principle, the gluons with transverse momentum Q s are extended 
in the longitudinal and proper time directions by ~ 1/Q S ] since the 
transverse area is ttR\, the four-volume is V4 ~ itR? a /Q 2 . The resulting 
four-density from (29) is then 

N 9 N g Q 2 4 

p ^ as v-r as ^~ Q - (31) 

where at the last stage we have used the non-linearity condition (30), 
Pi ~ Qs- It * s eas Y to see that (31) coincides with the saturation 
condition (26), since the number of gluons in the infinite momentum 
frame N g ~ xG(x, Q 2 ). 

In view of the significance of saturation criterion for the rest of the 
material in these lectures, let us present yet another argument, tradi- 
tionally followed in the discussion of classical limit in electrodynamics 
(Berestetskii,Lifshitz and Pitaevskii, 1982). The energy of the gluon 
field per unit volume is ~ E a2 . The number of elementary "oscillators 
of the field", also per unit volume, is ~ u; 3 . To get the number of the 
quanta in the field we have to divide the energy of the field by the 
product of the number of the oscillators ~ u; 3 and the average energy 
Huj of the gluon: 

The classical approximation holds when S> 1. Since the energy 
lo of the oscillators is related to the time At over which the average 
energy is computed by u ~ 1/At, we get 

^ » w?- (33) 

Note that the quantum mechanical uncertainty principle for the energy 
of the field reads 

E a2 uj 4 ~ h, (34) 
so the condition (33) indeed defines the quasi-classical limit. 
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Since E a2 is proportional to the action density p4, and the typ- 
ical time is At ~ using (31) we finally get that the classical 
description applies when 



kl < a s -^ = Ql (35) 



3.3. The absence of mini-jet correlations 

When the occupation numbers of the field become large, the matrix 
elements of the creation and annihilation operators of the gluon field 
defined by 

A" = y(cr At +ct At*) (36) 

ka ka ka ka' v ' 

k,a 

become very large, 

N ka = (ijkJ » 1- (37) 

so that one can neglect the unity on the r.h.s. of the commutation 
relation 

h ct - ct c t = 1 (38) 

ka ka ka ka V > 

and treat these operators as classical c— numbers. 

This observation, often used in condensed matter physics, espe- 
cially in the theoretical treatment of superfluidity, has important con- 
sequences for gluon production - in particular, it implies that the 
correlations among the gluons in the saturation region can be neglected: 

(A(k 1 )A(k 2 )...A(k n )) ~ (A(k 1 ))(A(k 2 ))...(A(k n )). (39) 

Thus, in contrast to the perturbative picture, where the produced mini- 
jets have strong back-to-back correlations, the gluons resulting from the 
decay of the classical saturated field are uncorrelated at k±_ 5; Q s . 

Note that the amplitude with the factorization property (39) is 
called point-like. However, the relation (39) cannot be exact if we con- 
sider the correlations of final-state hadrons - the gluon mini-jets can- 
not transform into hadrons independently. These correlations caused 
by color confinement however affect mainly hadrons with close three- 
momenta, as opposed to the perturbative correlations among mini-jets 
with the opposite three-momenta. 

It will be interesting to explore the consequences of the factorization 
property of the classical gluon field (39) for the HBT correlations of 
final-state hadrons. It is likely that the HBT radii in this case reflect 
the universal color correlations in the hadronization process. 
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Another interesting property of classical fields follows from the re- 
lation 

<(ct Cr f) - (cl Cr ) 2 = (ct Cr ), (40) 

which determines the fluctuations in the number of produced glu- 
ons. We will discuss the implications of Eq. (40) for the multiplicity 
fluctuations in heavy ion collisions later. 



4. Classical QCD in action 

4.1. Centrality dependence of hadron production 

In nuclear collisions, the saturation scale becomes a function of central- 
ity; a generic feature of the quasi-classical approach - the proportion- 
ality of the number of gluons to the inverse of the coupling constant 
(29) - thus leads to definite predictions (Kharzeev and Nardi, 2001) on 
the centrality dependence of multiplicity. 

Let us first present the argument on a qualitative level. At different 
centralities (determined by the impact parameter of the collision) , the 
average density of partons (in the transverse plane) participating in the 
collision is very different. This density p is proportional to the average 
length of nuclear material involved in the collision, which in turn ap- 
proximately scales with the power of the number N par t of participating 
1/3 

nucleons, p ~ N p ' art . The density of partons defines the value of the 
saturation scale, and so we expect 

Ql ~ K'arf (41) 

The gluon multiplicity is then, as we discussed above, is 

d V a s (QlY { > 

where Sa is the nuclear overlap area, determined by atomic number 
and the centrality of collision. Since Sa Q 2 s ~ Npart by definitions of 
the transverse density and area, from (42) we get 



dN, 



9 



dr) 



N paH ^ N part , (43) 



which shows that the gluon multiplicity shows a logarithmic deviation 
from the scaling in the number of participants. 

To quantify the argument, we need to explicitly evaluate the average 
density of partons at a given centrality. This can be done by using 
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Glauber theory, which allows to evaluate the differential cross section 
of the nucleus-nucleus interactions. The shape of the multiplicity dis- 
tribution at a given (pseudo) rapidity rj can then be readily obtained 
(see, e.g., (Kharzeev, Lourengo, Nardi and Satz, 1997)): 

= J d 2 bV{n-b) (l-Po(6)), (44) 

where Po(b) is the probability of no interaction among the nuclei at a 
given impact parameter b: 

P (b) = (l-a NN T AB (b)) AB ; (45) 

(Inn is the inelastic nucleon-nucleon cross section, and 1ab(&) is the 
nuclear overlap function for the collision of nuclei with atomic numbers 
A and B; we have used the three-parameter Woods-Saxon nuclear 
density distributions (De Jager, De Vries and De Vries, 1974). 
The correlation function V(n;b) is given by 

n 1 ( (n-n(b)) 2 \ 

T(n; b) = - exp — , 46 

here nib) is the mean multiplicity at a given impact parameter b; the 
formulae for the number of participants and the number of binary 
collisions can be found in (Kharzeev, Lourengo, Nardi and Satz, 1997). 
The parameter a describes the strength of fluctuations; for the classical 
gluon field, as follows from (40), a = 1. However, the strength of fluc- 
tuations can be changed by the subsequent evolution of the system and 
by hadronization process. Moreover, in a real experiment, the strength 
of fluctuations strongly depends on the acceptance. In describing the 
PHOBOS distribution (PHOBOS Coll., 2000), we have found that the 
value a = 0.6 fits the data well. 

In Fig. 4, we compare the resulting distributions for two different 
assumptions about the scaling of multiplicity with the number of par- 
ticipants to the PHOBOS experimental distribution, measured in the 
interval 3 < \r]\ < 4.5. One can see that almost independently of the- 
oretical assumptions about the dynamics of multiparticle production, 
the data are described quite well. At first this may seem surprising; the 
reason for this result is that at high energies, heavy nuclei are almost 
completely "black"; unitarity then implies that the shape of the cross 
section is determined almost entirely by the nuclear geometry. We can 
thus use experimental differential cross sections as a reliable handle on 
centrality. This gives us a possibility to compute the dependence of the 
saturation scale on centrality of the collision, and thus to predict the 
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Figure 4- Charged multiplicity distribution at \fs = 130 A GeV; solid line (his- 
togram) - PHOBOS result; dashed line - distribution corresponding to participant 
scaling (x = 0); dotted line - distribution corresponding to the 37% admixture of 
"hard" component in the multiplicity; see text for details. 

centrality dependence of particle multiplicities, shown in Fig. 5. (see 
(Kharzeev and Nardi, 2001) for details). 

4.2. Energy dependence 

Let us now turn to the discussion of energy dependence of hadron 
production. In semi-classical scenario, it is determined by the varia- 
tion of saturation scale Q s with Bjorken x = Q a /\fs. This variation, 
in turn, is determined by the x— dependence of the gluon structure 
function. In the saturation approach, the gluon distribution is related 
to the saturation scale by Eq.(26). A good description of HERA data is 
obtained with saturation scale Q 2 S = 1 4- 2 GeV 2 with W - dependence 
(W = y/s is the center-of-mass energy available in the photon-nucleon 
system) (Golec-Biernat and Wiisthof, 1999) 

Ql oc W\ (47) 

where A ~ 0.254-0.3. In spite of significant uncertainties in the determi- 
nation of the gluon structure functions, perhaps even more important 
is the observation (Golec-Biernat and Wiisthof, 1999) that the HERA 
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Figure 5. Centrality dependence of the charged multiplicity per participant pair 
near r\ = at y/s = 130 A GeV; the curves represent the prediction based on the 
conventional eikonal approach, while the diamonds correspond to the high density 
QCD prediction (see text). The square indicates the pp multiplicity. 



data exhibit scaling when plotted as a function of variable 

Qq \xo 



T = 



(48) 



where the value of A is again within the limits A ~ 0.25 0.3. In 
high density QCD, this scaling is a consequence of the existence of 
dimensionful scale (Gribov, Levin and Ryskin, 1983; McLerran and 
Venugopalan, 1994)) 

Q 2 s (x) = Ql (x /x) x . (49) 

Using the value of Q 2 S ~ 2.05 GeV 2 extracted (Kharzeev and Nardi, 
2001) at y/i = 130 GeV and A = 0.25 (Golec-Biernat and Wiisthof, 
1999) used in (Kharzeev and Levin, 2000), equation (59) leads to the 
following approximate formula for the energy dependence of charged 
multiplicity in central Au — Au collisions: 

/ 2 dN ^ ,0.87 ( ^ V ) r 5 v 



N, 



part 



T]<1 



3.93 + 0.25 In 



'yfl (GeV)~ 



130 



(50) 
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At y/s = 130 GeV, we estimate from Eq.(50) 2/N part dN cn /dr] | T) <i = 
3.42 ±0.15, to be compared to the average experimental value of 3.37 ± 
0.12 (PHOBOS Coll., 2000; PHENIX Coll., 2001; STAR Coll., 2001; 
BRAHMS Coll., 2001). At y/s = 200 GeV, one gets 3.91 ± 0.15, to be 
compared to the PHOBOS value (PHOBOS Coll., 2000) of 3.78 ±0.25. 
Finally, at y/s = 56 GeV, we find 2.62 ± 0.15, to be compared to 
(PHOBOS Coll., 2000) 2.47±0.25. It is interesting to note that formula 
(50), when extrapolated to very high energies, predicts for the LHC 
energy a value substantially smaller than found in other approaches: 



dN, 



ch 



N, 



part 



df] 



10.8 ± 0.5; = 5500 GeV, 



(51) 



7)<1 



corresponding only to a factor of 2.8 increase in multiplicity between 
the RHIC energy of y/s = 200 GeV and the LHC energy of y/s = 
5500 GeV (numerical calculations show that when normalized to the 
number of participants, the multiplicity in central Au — Au and Pb — Pb 
systems is almost identical) . The energy dependence of charged hadron 
multiplicity per participant pair is shown in Fig. 6. 



Multiplicity per participant pair 

C.m.s. energy, GeV 




50 100 150 200 250 

Figure 6. Energy dependence of charged multiplicity per participant pair at RHIC 
energies; solid line is the result (50). 
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One can also try to extract the value of the exponent A from the 
energy dependence of hadron multiplicity measured by PHOBOS at 
yfs = 130 GeV and at at = 56 GeV; this procedure yields A ~ 0.37, 
which is larger than the value inferred from the HERA data (and is very 
close to the value A ~ 0.38, resulting from the final-state saturation 
calculations (Eskola, Kajantie, Ruuskanen and Tuominen, 2000)). 

4.3. Radiating the classical glue 

Let us now proceed to the quantitative calculation of the (pseudo-) 
rapidity and centrality dependences (Kharzeev and Levin, 2001). We 
need to evaluate the leading tree diagram describing emission of gluons 
on the classical level, see Fig. 7 5 . 




Figure 7. The Mueller diagram for the classical gluon radiation. 

Let us introduce the unintegrated gluon distribution ipj± (x, fcf ) which 
describes the probability to find a gluon with a given x and transverse 
momentum kf inside the nucleus A. As follows from this definition, the 

a Note that this "mono-jet" production diagram makes obvious the absence of 
azimuthal correlations in the saturation regime discussed above, see eq (39). 
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unintegrated distribution is related to the gluon structure function by 

/Pt 
dk 2 ip A (x,k 2 ); (52) 

when p\ > Q 2 , the unintegrated distribution corresponding to the 
bremsstrahlung radiation spectrum is 

<p A (x,tf) ~ — i. (53) 

7T kf 

In the saturation region, the gluon structure function is given by (27); 
the corresponding unintegrated gluon distribution has only logarithmic 
dependence on the transverse momentum: 

<Pa(x, kf) ~ — ; k 2 <Q 2 , (54) 
a s 

where Sa is the nuclear overlap area, determined by the atomic numbers 
of the colliding nuclei and by centrality of the collision. 

The differential cross section of gluon production in a AA collision 
can now be written down as (Gribov, Levin and Ryskin, 1983; Gyulassy 
and McLerran, 1997) 

E ~7T = ~2 / dk t a s <PA(xi,kt) pa(x2, (p ~ k)t), (55) 



d 3 p N 2 -l p 2 

where x\^ = {pt/ \fs) exp(±r/), with r] the (pseudo)rapidity of the pro- 
duced gluon; the running coupling a s has to be evaluated at the scale 
Q 2 = max{k 2 , (p — k) 2 }. The rapidity density is then evaluated from 
(55) according to 

dN 1 f , 9 / „ da \ 



dy oaa 



where gaa is the inelastic cross section of nucleus-nucleus interaction. 

Since the rapidity y and Bjorken variable are related by In 1/x = y, 
the x— dependence of the gluon structure function translates into the 
following dependence of the saturation scale Q 2 on rapidity: 

Q 2 (s; ±y) = Q 2 s (s; y = 0) exp(±Ay). (57) 

As it follows from (57), the increase of rapidity at a fixed W = 
y/s moves the wave function of one of the colliding nuclei deeper into 
the saturation region, while leading to a smaller gluon density in the 
other, which as a result can be pushed out of the saturation domain. 
Therefore, depending on the value of rapidity, the integration over the 
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transverse momentum in Eqs. (55), (56) can be split in two regions: 
i) the region Aqcd < kt < Qs,min m which the wave functions are 
both in the saturation domain; and ii) the region A << Q S) min < 
k t < Qs. max in which the wave function of one of the nuclei is in the 
saturation region and the other one is not. Of course, there is also the 
region of k t > Q s , ma x, which is governed by the usual perturbative 
dynamics, but our assumption here is that the role of these genuine 
hard processes in the bulk of gluon production is relatively small; in 
the saturation scenario, these processes represent quantum fluctuations 
above the classical background. It is worth commenting that in the 
conventional mini-jet picture, this classical background is absent, and 
the multi-particle production is dominated by perturbative processes. 
This is the main physical difference between the two approaches; for 
the production of particles with pt >> Q s they lead to identical results. 

To perform the calculation according to (56), (55) away from y = 
we need also to specify the behavior of the gluon structure function 
at large Bjorken x (and out of the saturation region). At x — > 1, this 
behavior is governed by the QCD counting rules, xG(x) ~ (1 — x) 4 , so 
we adopt the following conventional form: xG(x) ~ x~ x (1 — x) 4 . 

We now have everything at hand to perform the integration over 
transverse momentum in (56), (55); the result is the following (Kharzeev 
and Levin, 2001): 

^ = const S A Q\ mm In (||^) x 

. + 2 lQ {Ql min ) V yfi ) 

where the constant is energy-independent, Sa is the nuclear overlap 
area, Q 2 S = Q^(s;y = 0), and Q s ^ m in(max) are defined as the smaller 
(larger) values of (57); at y = 0, Q\ min = Q% max = QsOO = Qs( s o) x 
x (s/so) A / 2 . The first term in the brackets in (58) originates from the re- 
gion in which both nuclear wave functions are in the saturation regime; 
this corresponds to the familiar ~ (l/«s) Q 2 S R\ term in the gluon 
multiplicity. The second term comes from the region in which only one 
of the wave functions is in the saturation region. The coefficient 1/2 in 
front of the second term in square brackets comes from kt ordering of 
gluon momenta in evaluation of the integral of Eq.(55). 

The formula (58) has been derived using the form (54) for the 
unintegrated gluon distributions. We have checked numerically that 
the use of more sophisticated functional form of (fA taken from the 
saturation model of Golec-Biernat and Wusthoff (Golec-Biernat and 



(58) 
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Wusthof, 1999) in Eq.(55) affects the results only at the level of about 
3%. 

Since SaQ 2 s ~ N part (recall that Q 2 S » ^qcd ^ s defined as the 
density of partons in the transverse plane, which is proportional to the 
density of participants), we can re-write (58) in the following final form 
(Kharzeev and Levin, 2001) 

(*£)- A|! " x 



dN 
dy 



c N. 



part 



±V P -Mv\ 
so 



1 + A|y| 1 



Q± e (l+A/2)|y| 



(59) 



with Qg(s) = Q1(sq) (s/sq) x / 2 . This formula is the central result of our 
paper; it expresses the predictions of high density QCD for the energy, 
centrality, rapidity, and atomic number dependences of hadron multi- 
plicities in nuclear collisions in terms of a single scaling function. Once 
the energy-independent constant c ~ 1 and Q 2 (so) are determined at 
some energy so, Eq. (59) contains no free parameters. At y = the 
expression (58) coincides exactly with the one derived in (Kharzeev 
and Nardi, 2001), and extends it to describe the rapidity and energy 
dependences. 



4.4. Converting gluons into hadrons 



The distribution (59) refers to the radiated gluons, while what is mea- 
sured in experiment is, of course, the distribution of final hadrons. We 
thus have to make an assumption about the transformation of gluons 
into hadrons. The gluon mini-jets are produced with a certain virtu- 
ally, which changes as the system evolves; the distribution in rapidity 
is thus not preserved. However, in the analysis of jet structure it has 
been found that the angle of the produced gluon is remembered by 
the resulting hadrons; this property of "local parton-hadron duality" 
(see (Dokshitzer, 1998) and references therein) is natural if one as- 
sumes that the hadronization is a soft process which cannot change 
the direction of the emitted radiation. Instead of the distribution in 
the angle 6, it is more convenient to use the distribution in pseudo- 
rapidity r] = — lntan(0/2). Therefore, before we can compare (58) to 
the data, we have to convert the rapidity distribution (59) into the gluon 
distribution in pseudo-rapidity. We will then assume that the gluon 
and hadron distributions are dual to each other in the pseudo-rapidity 
space. 
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To take account of the difference between rapidity y and the mea- 
sured pseudo-rapidity 77, we have to multiply (58) by the Jacobian of 
the y <-> r\ transformation; a simple calculation yields 



h(ri;pt;m) 



cosh 7/ 



m 2 +pi 



(60) 



+ sinh 2 r\ 



where m is the typical mass of the produced particle, and pt is its 
typical transverse momentum. Of course, to plot the distribution (59) 
as a function of pseudo-rapidity, one also has to express rapidity y in 
terms of pseudo-rapidity 77; this relation is given by 



y{v,Pf,m) = - In 



+ sinh i] + sinh rj 



1 t Pt + sinh 2 77 — sinh t? 
PI . 



(61) 



obviously, h(r];p t ;m) = dy(rj;pt;m)/drj. 

We now have to make an assumption about the typical invariant 
mass 777 of the gluon mini-jet. Let us estimate it by assuming that the 
slowest hadron in the mini-jet decay is the p-resonance, with energy 
E p = (to 2 +p 2 p j +P 2 p ,z) 1 ^ 2 1 where the z axis is pointing along the mini- 
jet momentum. Let us also denote by Xj the fractions of the gluon 
energy q$ carried by other, fast, i particles in the mini-jet decay. Since 
the sum of transverse (with respect to the mini-jet axis) momenta of 
mini-jet decay products is equal to zero, the mini-jet invariant mass to 
is given by 



m )et 



777 



~ 2^2xiq z -(mp :t -pp iZ ) = 2Q s -m ef f, (62) 

i 

where 777 Pit = (m 2 +p 2 it ) 1 ^ 2 - In Eq. (62) we used that J2i x i = 1 an d 
Qo ~ Qz = Qs- Taking p p<z « p Pjt w 300 MeV and p mass, we obtain 
777 e// w 0.5 GeV. 

We thus use the mass m 2 ~ 2Q s m e ff ~ Q s ■ 1 GeV in Eqs. (60,61). 
Since the typical transverse momentum of the produced gluon mini-jet 
is Q s , we take pt = Q s in (60). The effect of the transformation from 
rapidity to pseudo-rapidity is the decrease of multiplicity at small 77 
by about 25 — 30%, leading to the appearance of the ^ 10% dip in the 
pseudo-rapidity distribution in the vicinity of 7/ = 0. We have checked 
that the change in the value of the mini-jet mass by two times affects 
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Figure 8. Centrality dependence of charged hadron production per participant at 
different pseudo-rapidity r\ intervals in Au — Au collisions at i/s = 130 GeV; from 
(Kharzeev and Levin, 2001), the data are from (PHOBOS Coll., 2000). 



the Jacobian at central pseudo-rapidity to about ~ 10%, leading to 
~ 3% effect on the final result. 

The results for the Au—Au collisions at y/s = 130 GeV are presented 
in Figs 8 and 9. In the calculation, we use the results on the dependence 
of saturation scale on the mean number of participants at y/s = 130 
GeV from (Kharzeev and Nardi, 2001), see Table 2 of that paper. The 
mean number of participants in a given centrality cut is taken from 
the PHOBOS paper (PHOBOS Coll., 2000). One can see that both the 
centrality dependence and the rapidity dependence of the y/s = 130 
GeV PHOBOS data are well reproduced below r] ~ ±4. The rapidity 
dependence has been evaluated with A = 0.25, which is within the 
range A = 0.25 -j- 0.3 inferred from the HERA data (Golec-Biernat 
and Wusthof, 1999). The discrepancy above rj ~ ±4 is not surprising 



2( dN ch /dr|)/N part 
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Figure 9. Pseudo-rapidity dependence of charged hadron production at different 
cuts on centrality in Au — Au collisions at i/s = f30 GeV; from (Kharzeev and 
Levin, 2001), the data are from (PHOBOS Coll., 2000). 



since our approach does not properly take into account multi-parton 
correlations which are important in the fragmentation region. 

Our predictions for Au — Au collisions at y/s = 200 GeV are pre- 
sented in (Kharzeev and Levin, 2001). The only parameter which gov- 
erns the energy dependence is the exponent A, which we assume to be 
A ~ 0.25 as inferred from the HERA data. The absolute prediction 
for the multiplicity, as explained above, bears some uncertainty, but 
there is a definite feature of our scenario which is distinct from other 
approaches. It is the dependence of multiplicity on centrality, which 
around rj = is determined solely by the running of the QCD strong 
coupling (Kharzeev and Nardi, 2001). As a result, the centrality depen- 
dence at \/s = 200 GeV is somewhat less steep than at -y/i = 130. While 
the difference in the shape at these two energies is quite small, in the 
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perturbative mini-jet picture this slope should increase, reflecting the 
growth of the mini-jet cross section with energy (Wang and Gyulassy, 
2001). 

4.5. Further tests 

Checking the predictions of the semi-classical approach for the cen- 
trality and pseudo-rapidity dependence at \/s = 200 GeV is clearly 
very important. What other tests of this picture can one devise? The 
main feature of the classical emission is that it is coherent up to the 
transverse momenta of about \/2 Q s (about ~ 2 GeV/c for central 
Au — Au collisions). This means that if we look at the centrality depen- 
dence of particle multiplicities above a certain value of the transverse 
momentum, say, above 1 GeV/c, it should be very similar to the depen- 
dence without the transverse momentum cut-off. On the other hand, in 
the two-component "soft plus hard" model the cut on the transverse 
momentum would strongly enhance the contribution of hard mini-jet 
production processes, since soft production mechanisms presumably 
do not contribute to particle production at high transverse momenta. 
Of course, at sufficiently large value of the cutoff all of the observed 
particles will originate from genuine hard processes, and the centrality 
dependence will become steeper, reflecting the scaling with the number 
of collisions. It will be very interesting to explore the transition to this 
hard scattering regime experimentally. 

Another test, already discussed above (see eq.(39)) is the study of 
azimuthal correlations between the produced high p t particles. In the 
saturation scenario these correlations should be very small below p t ~ 2 
GeV/c in central collisions. At higher transverse momenta, and/or for 
more peripheral collisions (where the saturation scale is smaller) these 
correlations should be much stronger. 



5. Does the vacuum melt? 

The approach described above allows us to estimate the initial energy 
density of partons achieved at RHIC. Indeed, in this approach the 
formation time of partons is To — 1 /Q s , and the transverse momenta of 
partons are about kt — Q s . We thus can use the Bjorken formula and 
the set of parameters deduced above to estimate (Kharzeev and Nardi, 
2001) 
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Figure 10. Distributions of instanton sizes in vacuum for QCD with three light 
flavors (upper curve) versus the distribution of instanton sizes in the saturation 
environment produced by a collision of two identical nuclei for c = 1 (middle curve) 
and c = 2 In 2 (lower curve) with Q 2 a = 2GeV 2 ; from (Kharzeev, Kovchegov and 
Levin, 2002). 



for central Au — Au collisions at ^/s = 130 GeV. This value is well 
above the energy density needed to induce the QCD phase transition 
according to the lattice calculations. However, the picture of gluon 
production considered above seems to imply that the gluons simply flow 
from the initial state of the incident nuclei to the final state, where they 
fragment into hadrons, with nothing spectacular happening on the way. 
In fact, one may even wonder if the presence of these gluons modifies 
at all the structure of the physical QCD vacuum. 

To answer this question theoretically, we have to possess some knowl- 
edge about the non-perturbative vacuum properties. While in general 
the problem of vacuum structure still has not been solved (and this 
is one of the main reasons for the heavy ion research!), we do know 
one class of vacuum solutions - the instantons. It is thus interesting to 
investigate what happens to the QCD vacuum in the presence of strong 
external classical fields using the example of instantons (Kharzeev, 
Kovchegov and Levin, 2002). 

The problem of small instantons in a slowly varying background 
field was first addressed in (Callan, Dashen and Gross, 1979; Shifman, 
Vainshtein and Zakharov, 1980) by introducing the effective instanton 
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lagrangian L^j{x) 

L lff( x o) = J dpn (p) dR exp {-—- p 2 rj™ v R aa ' G< (x )) (64) 

in which no(p) is the instanton size distribution function in the vacuum, 
VaLu is t ne 't Hooft symbol in Minkowski space, and R aa is the matrix 
of rotations in color space, with dR denoting the averaging over the 
instanton color orientations. 

The complete field of a single instanton solution could be recon- 
structed by perturbatively resumming the powers of the effective in- 
stanton lagrangian which corresponds to perturbation theory in powers 
of the instanton size parameter p 2 . In our case here the background 
field arises due to the strong source current J". The current can be 
due to a single nucleus, or resulting from the two colliding nuclei. Per- 
turbativc resummation of powers of the source current term translates 
itself into resummation of the powers of the classical field parameter 
a 2 A 1 / 3 (McLerran and Venugopalan, 1994; Kovchegov, 1996). Thus 
the problem of instantons in the background classical gluon field is 
described by the effective action in Minkowski space 

Seff = J d 4 x (-\d%(x)G^(x) + Li ff (x) + lJ eff (x) + J°^(z)) . 

(65) 

The problem thus is clearly formulated; by using an explicit form 
for the radiated classical gluon field, it was possible to demonstrate 
(Kharzeev, Kovchegov and Levin, 2002) that the distribution of in- 
stantons gets modified from the original vacuum one no(p) to 

n£W = ».M exp (- 8g // ( | To)2 ) , (66) 

where To is the proper time. The result Eq. (66) shows that large size 
instantons are suppressed by the strong classical fields generated in the 
nuclear collision (see Fig. 10) 6 . The vacuum does melt! 

The results presented here were obtained together with Hirotsugu 
Fujii, Yuri Kovchegov, Eugene Levin, and Marzia Nardi. I am very 
grateful to them for the most enjoyable collaboration. I also wish to 
thank Jean-Paul Blaizot, Yuri Dokshitzer, Larry McLerran, Al Mueller 
and Raju Venugopalan for numerous illuminating discussions on the 

6 Of course, at large proper times To — » oo the vacuum "cools off", and the 
instanton distribution returns to the vacuum one. 
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